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1 Introduction 

Q^ ' In this paper we study a class of nonlinear Schrodinger equations which admit families of 

small solitary wave solutions. We consider solutions which are small in the energy space 
H , and decompose them into solitary wave and dispersive wave components. The goal is 
to establish the asymptotic stability of the solitary wave and the asymptotic completeness 
of the dispersive wave. That is, we show that as t — > oo, the solitary wave component 
converges to a fixed solitary wave, and the dispersive component converges to a solution of 
the free Schrodinger equation. 

Let us briefly supply some background. Solutions of dispersive partial differential equa- 
tions (with repulsive nonlinearities) tend to spread out in space, although they often have 
conserved L 2 mass. There has been extensive study of this phenomenon, usually referred to 
as scattering theory. These equations include Schrodinger equations, wave equations, and 
KdV equations. However, such equations can also possess solitary wave solutions which have 
localized spatial profiles that are constant in time (e.g., if the nonlinearity is attractive or if 
a linear potential is present). To understand the asymptotic dynamics of general solutions, 
it is essential to study the interaction between the solitary waves and the dispersive waves. 
The matter becomes more involved when the linearized operator around the solitary wave 
possesses multiple eigenvalues, which correspond to excited states. The interaction between 
eigenstates (mediated by the nonlinearity) is very delicate, and few results are known. 
For nonlinear Schrodinger equations with solitary waves, there are three types of results: 

1. Control of solutions in a finite time interval, and construction of all-time solutions 
with specified asymptotic behaviors {scattering solutions, see I10j). This type of result 
does not allow sufficient time for the excited state interaction to make a difference. 

2. Orbital stability of solitary waves. A solution stays close to the family of nonlinear 
bound states if it is initially close. This is usually proved by energy arguments, see e.g. 

123 III EEH1- 

3. Asymptotic stability of solitary waves. Here, one must assume that the spectrum 
of the linearized operator enjoys certain spectral properties (for example, has only one 
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eigenvalue, or has multiple "well-placed" eigenvalues). Furthermore, the initial data are 
typically assumed to be localized, so that the dispersive wave has fast local decay. Even 
under restrictive spectral assumptions, only perturbation problems can be treated for large 
solitary waves (see [HHH, also HJ IS] for 1-D results), while more general results can be 
obtained for small solitary waves |27I l lT7 | l2l^l2^ l231 l2fi l l27 | l23]. 
In this paper, we study small solutions of the equation 

idti/) = (-A + V)il> + gty), iP(0,-)=tP £H 1 (M 3 ) (1.1) 

with small data: ||"0q|| if i << 1 (this is equivalent to considering a nonlinearity multiplied 
by a small constant). Although we only consider the problem for x G M 3 , the results and 
methods can be extended to spatial dimensions d > 3. 

Here, g(ip) is either a pointwise nonlinearity or a Hartree-type (non-local) nonlinearity 
(or their sum), satisfying gauge covariance: 

g(ije ia )=g^)e ia , with g{\^\) G R. (1.2) 

More detailed assumptions are given below. In either case, we can find a functional G : 
H 1 — > M, satisfying G{ifje %a ) = G(ip) (gauge invariance), and 

d° £ G(*P + erf) := j-G{^ + e V )\ e=0 = Re^), r,). (1.3) 

Here we denote the inner product in L 2 by 

(a,b) := / abdx. (1.4) 

Jr3 

Under suitable assumptions, the L 2 -norm HVKOIIl 2 an d the Hamiltonian 

I [ (\ViP\ 2 + V\iP\ 2 )dx + GW (1.5) 

are constant in time. Using these conserved quantities, and the smallness of HV'oIIhi; one 
can prove a uniform estimate sup t IIV'WIIhi ^ 1 an d obtain global well-posedness. This is, 
however, not used in our proof. 

We do not assume that is localized (ipo G L 1 (M 3 ), for example, or in a weighted 
space), as is usually done. As a result, we cannot expect a time decay rate for L p -norms 
(p > 2) of the dispersive part of the solution. However, the space H 1 is natural, as it is 
intimately related to the Hamiltonian structure, and persists globally in time (in contrast 
to weighted spaces, whose smallness persists only for short time due to dispersion, and L 1 , 
which may be instantaneously lost and so does not seem to have physical relevance). A 
related motivation comes from the situation where the linearized operator around a solitary 
wave has many eigenvalues. In this case, the dispersive component tends to decay very 
slowly. It is thus essential to be able to remove the localization assumption on the data. 
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We assume that — A + V supports only one eigenvalue eo < 0, which is non-degenerate, 
and we denote by (fio the corresponding positive, normalized eigenfunction. More detailed 
assumptions on V are given below. Under these assumptions, there exists a family of small 
"nonlinear bound states" Q = Q[z], parameterized by small z = (4>o,Q) € C, which satisfy 
Q[z] — z(j)Q = o{z) _L </>o, and solve the nonlinear eigenvalue problem 

(-A + V)Q + g(Q) = EQ, E = E[z] = e + o(z) G R. (1.6) 

See Lemma 12. II for details. Gauge covariance is inherited by Q: 

Q[ze ia ] = Q[z]e ia , (1.7) 

and so E[z] = E[\z\]. The nonlinear bound states give rise to exact solitary-wave solutions 
ip(x,i) = Q(x)e~ tEt of (|1.1|) . Q[z] is differentiable in z if we regard it as a real vector 

z = zi + iz 2 <-» (zi,z 2 ) € R 2 . (1.8) 

We will denote its z-derivatives by 

DxQ[z] := -^Q[z], D 2 Q[z] := —Q[z] (1.9) 

(we use the symbol D in order to distinguish them from space or time derivatives). Then 
DQ[z] denotes the Jacobian matrix, regarded as a R- linear map on C: 

£>Q[z]:C->C, DQ[z]w ^ D 1 Q[z]Rew + iD 2 Q[z]lmw. (1.10) 

The gauge covariance of Q[z] implies that 

DQ[z}iz = iQ[z]. (1.11) 

Given a general solution ?/>(£) of ()1.1|) . it is natural to decompose it into solitary wave 
and dispersive wave components: 

m = Qim]+v(t). (i.i2) 

For any such decomposition, (|1.1|) yields an equation for ij: 

if] = H[z]r] + E[z]Q[z] - iDQ[z]z + F 2 (z, rj), (1.13) 
where i7[z] denotes the linearized operator around Q[z], 

H[z]ri := (-A + V) V + d° £ g(Q + er/), (1.14) 
and F 2 collects terms which are higher-order in rj: 

F 2 (z, V ) := g(Q + v ) - g(Q) - d° £ g(Q + erj). (1.15) 
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The decomposition Q1.12|) is of course not unique. To specify the path z(t) uniquely, 
we impose an orthogonality condition which will make r/ dispersive. Since the linearization 
destroys gauge invariance, the linearized operator H[z] is not complex-linear. It is, however, 
symmetric if we regard C as M 2 , and use the reduced inner product: 

(a, b) := Re(a, b) = / (Re a Re 6 + Imalmfr) dx. (1-16) 

Jm.3 

The symmetry of H[z] follows from Q1.14[) and 

(d° e g(Q + £77), C> = d° e (g(Q + er,), () = d° £ d° s G(Q + er? + SQ- (1.17) 

We define the "continuous spectral subspace" 

K c [z] := {i)£L 2 : (i V , D l Q[z\) = D 2 Q[z}) = 0} . (1.18) 

This is an invariant subspace of i(H[z] — E[z]), as follows from the relation 

(H[z\ - E^DjQlz] = iPtEWQlz] (1.19) 

(which is the result of differentiating (|1.6|) ). together with (|1.11|) . Restricting to "K c [z\ 
eliminates non-decaying solutions of the linear equation c\r/ = —i(H[z] — E[z])r) for fixed z. 
When zeK+K c [z] is just the orthogonal complement of {Q^^^Q} m tb e inner product 
(■, •} (this subspace is often encountered in the literature). 

As we will show in Lemma l'2. 31 we can uniquely decompose ip(t) as 

rP(t) = Q[z(t)} + r?(t), r,(t) G H c [z(t)}. (1.20) 

The requirement rj(t) E ^K c [z(t)] determines z(t) uniquely. An evolution equation for z(t) 
is derived from differentiating the relation (irj, DjQ[z]) = with respect to t, and using 
equation 1)1. 13|) (see (|3.12|) 'l. Our goal is to prove the asymptotic stability of Q[z(t)] and 
the asymptotic completeness of r)(t). 

We now state precise assumptions on the potential V, and on the nonlinearity g. We 
denote the usual Lorentz space by L p,q = (L°°, L 1 )i/ p q for 1 < p < oo and 1 < q < oo (see 
0). W 1,p denotes the usual Sobolev space. 

Assumption 1: V is a real- valued function belonging to L 2 + L°° . (We note that under 
this assumption, —A+V is a self-adjoint operator on L 2 , with domain H 2 . See, eg. |18|). Its 
negative part V- := max{0, — V} is further assumed to satisfy 1 1 1 1 ^ .j.^ 00 ) (-{ la; I >i2}) ~^ ® as 
R — > oo. We suppose — A + V has only one eigenvalue eo < 0, and let </>o be a corresponding 
normalized eigenvector, eo is simple and cfio can be taken to be positive (,18;). Denote the 
projections onto the discrete and continuous spectral subspaces of —A + V by 

P d = M<t>ar), P c = l-P d . (1.21) 
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The following Strichartz estimates are assumed to hold: 



e it(A-V)p c , 
* Js(A-V) 



X 



< 



H 1 : 



P c F(s)ds 



< \\F\ 



(1.22) 



L a w i,B/5 , 



X 



where X := LfH 1 n L 2 W 1,6 n L 2 L 6 ' 2 . 

We remark that the Strichartz estimates of Assumption 1 hold when, e.g., 



|V»| < C(l + \x\ 



-3-e 



11-23) 



for some e > 0, and the bottom of the continuous spectrum, zero, is neither an eigenvalue 
nor a resonance. Estimates without derivatives can be proved by applying the L 1 -L°° decay 
estimate EDI EES] to the endpoint Strichartz estimate ^H] , where the stronger estimate 
in the Lorentz space was actually proved. We need the Lorentz space L 6 ' 2 estimate only to 
handle the critical case of the Hartree equation (with convolution potentials decaying like 
l/|x| 2 ). Estimates of the derivatives can be obtained by using the equivalence 



I LP 



H 



-a + v + iiyii +i 



(1.24) 



for 1 < p < oo, and the commutativity of e lt ^~ v ^ with H 1 / 2 . The equivalence can be shown 
by applying the complex interpolation for fractional powers [221 §1.15.3] to the equivalence 
in W 2 ' p ', using the boundedness of imaginary powers H ls , which follows by [Jj from the fact 
that e~ tH is a positivity preserving contraction semigroup on L p . 

Assumption 2: The nonlinearity is assumed to be of one of the following two forms, 
or their sum: 

(a) g : C — > C is a function satisfying gauge covariance (J1.2|) which, when restricted to 
M, is twice differentiable, with g(0) = g'(0) = 0, and 



(b) g(ip) = (<& * l^l 2 )^, where $ is a real potential, and 

$ G L 1 + L 3 / 2 >°°. 



(1.25) 



(1.26) 



Examples of nonlinearities satisfying Assumption 2 include 

d s 
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(i/>) =a|V| 4/3 ^ + 6|V|V+ [(] 



1 3-e 



+ 



(1.27) 



where a, 6, c, d G R, < e < 1. 

We can now state our main theorems. 
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Theorem 1.1 (Asymptotic stability and completeness) Let Assumptions 1 and 2 hold. 
Every solution ip of with data ipQ sufficiently small in H 1 can be uniquely decomposed 
as 

i>{t)=Q[z(t)]+r ] {t), (1-28) 
with differentiable z(t) € C and rj(t) £ !K c [z(t)] satisfying 

WvWtfw^nLFH 1 + \\ z \\l?° ~ llV'oll/fi, 

(1.29) 

\\z + iE[z]z\\ L i Lr < Uo\\ 2 H i . 



Moreover, there exist > with — |^(0) | < HV'oll^i; and i] + G H 1 n RanP c such 
that 

W)\ moo, r?(t) - e ft ( A - v )»7 + , - (1.30) 



as t — > oo. 



The corresponding result with no bound state was obtained in |21II15| for small H 1 data 
and ^21 f° r large data with no potential and g(ip) = +\ip\ m ~ 1 ip. Results similar to Theorem 
11.11 in the case of localized initial data ifio and g(ip) = Xl'ipl 171 ^ 1 '^ were first obtained for the 
case ||f?(0)||^i nL i <C 1-2(0)1 by Soffer and Weinstein [20], and then extended to all tpo small 
in H 1 and weighted L 2 -spaces by Pillet and Wayne |17| . The latter work was extended to 
the 1-D case in j2H]. In all [201 El EEj- the solutions ij)(t) are decomposed with respect to 
fixed self-adjoint linear operators. A time-dependent decomposition similar to ()1.20[) seems 
to have first appeared in [3]. 

Theorem 1.2 (Nonlinear wave operator) Let Assumptions 1 and 2 hold. There exists 
6 > such that for any £ [0,6] and r/ + € H 1 n RanP c with \\t]+\\ h i < 6, there is a 
global solution ip(t) of (|l.ljl satisfying the conclusion of Theorem li.il with the prescribed 
asymptotic data and r] + . 

Special cases of Theorem 11.21 further assuming 3> ||7/+||#i n £i or moo = with 
||?7_)_|| f/ -i nL i -C 1 for g(ip) = ±\ip\ 2 ip, were obtained in |2H 

In Theorems II .11 and II .2\ one may replace e lt ^~ v ^ri + by e lt ^rj+ with rj+ € H 1 if asymp- 
totic completeness in H 1 of the wave operator between — A + V and —A holds. It holds, 



for example, if ([1.23ft holds. 

Theorem 11.11 implies in particular that any small solution looks like a solitary wave for 
large time locally in space. But due to the fact that the data is not assumed to be localized, 
we cannot, in general, have a convergence rate for it. In fact we have 

Theorem 1.3 (Examples of slow decay of dispersion) Let Assumptions 1 and 2 hold. 
For any nonempty ball 8cl 3 , there exists 6 > for which the following holds. For any 
positive function f(t) which goes to zero as t — > oo ; and any G [0,6], there exists a 
solution ip(t) of (jl.lj) satisfying the conclusions of Theorem M.ll and 

limsup ( inf \\m - Q[z']\\ mm /f(t)) = oo. (1.31) 

t^co \\Z |<2<5 > J 
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Remarks: Note that we impose a time-dependent condition rj(t) G J£ c [z(i)] instead of 
simpler conditions such as rj{t) G 5C C [0], i-e., (rj(t),(/>o) = (which is used in ^1|2H])- The 
reason is the following. If we assume (ji(t), 4>q) = 0, then the equations for z + iEz yield 

\z + iEz\ < 1(00,^)1 + 1(00,^)1, 

where A is some linear operator. The term (0q, Arj) is linear in r\ and hence is not integrable 
in time, in light of the estimate r\ 6 L^W 1 '®. Thus we cannot conclude that \z\ and E[z] have 
limits as t — > oo. This term drops out if we require rj(t) G 5C C [ Z ], an d the equation for z+iEz 
(and hence ^\z\) becomes quadratic in rj. Even without this term, we are forced to use an 
Lj-type Strichartz estimate in order to get convergence of \z\, since we cannot have better 
decay as long as we start with H 1 initial data. Thus the end point Strichartz estimates 
(|1.22j) are irreplaceable in our argument. As a bonus, they allow us to treat borderline 
nonlinearities such as g(tjj) = ±|^| 4 V> and g(ijj) = ±|^| 4 / 3 ^, which are not covered in the 
previous works |2L)l I17j . By considering the problem in higher Sobolev spaces H s , s > 1 , 
one may treat higher power nonlinearities g(ip) = zbl^l" 1 " 1 ^ for m G (5, 1 + 4/(3 — 2s)). 



2 Preliminaries 

In this section, we give three lemmas concerning the nonlinear ground states. 

Lemma 2.1 (Nonlinear ground states) There exists 5 > such that for each z G C 
with \z\ < 5, there is a solution Q[z] G H 2 n W 1 ' 1 of (jl.6|) with E = E[\z\] G R such that 

Q[z] = Z0 O + q[z], (q,(j) ) = 0. 

The pair (q, zE) is unique in the class 

\\q\\h 2 ^ ^ \E ~ e ol < 
Moreover, Q[ze' ia ] = Q[z]e ta , Q[\z\\ is real, and 

q[z] = o(z 2 ) ) oil 

DQ[z] = (l,i)<l> + o(z), D 2 Q[z] = o(l) J (2.1) 
E[z]=e + o(z), DE[z]=o(l), 

as z — > 0. 

A special case of this lemma is proved in |2()j , referring to . We will prove the lemma 
under weaker assumptions on V and g in the Appendix. 

The following is an immediate but useful corollary of this lemma. 

Lemma 2.2 (Continuous spectral subspace comparison) There exists 5 > such 
that, for each z G C with \z\ < 5, there is a bijective operator R[z] : IK c [0] — > "K c [z\ 
satisfying 

PcU [z] = R[z}- 1 . (2.2) 
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Moreover, R[z] — I is compact and continuous in z in the operator norm on any space Y 
satisfying H 2 n W 1 * 1 C Y C H~ 2 + L°° . 

We remark that no corresponding statement holds for the case of large solitary waves. 
Proof. R[z] is given by 

R[z]rj = 7] + (p a[z]ri, (2.3) 
where the operator a[z] : !K C [0] — ► C is defined by solving the equations 

{ V + <t> a[z]ri,D j Q[z})=Q, j = 1,2. (2.4) 

This is solvable due to the property DjQ[z] = (l,i)<po + o(z). Then R[z] is obviously the 
inverse of P c restricted onto !K c [z]. For any rj £ !K C [0] we have 

Hz] V \ < \(v,DQ[z})\ < o(z)\\ti\\h-'+l°°, (2-5) 

which implies compactness of R{z] — I in Y. The continuity in z follows from that of DQ[z\. 
□ 

Lemma 2.3 (Best decomposition) There exists S > such that any ip G H 1 satisfying 
\ H i < 8 can be uniquely decomposed as 



1> = Q[z]+rt, (2.6) 

where z S C, n G "K c [z] and \z\ + \\r]\\ H i < HV'lli? 1 - 

Proof. We look for a unique solution z of the equation A(z) = 0, where we define 

Aj{z) := (i{i> - Q[z]), DjQ[z]), j = 1, 2. 

Let 

n -~ llV'llfl'i- The Jacobian matrix of the map z *— > ^4(^) is written as 



DjA k (z) = {-iD 3 Q[z],D k Q[z\) + (ity - Q[z]), D^fcQfz]) 
= j — k + o(n + \z\), 



(2.7) 



by Lemma 12. II Let zq := (<po,ip). So \zq\ < n. Then from Lemma 12. II we have 

A(zq) = {i{i> - zqcPq) + o(n 2 ), (1, i)<t> Q + o(n)) = o(n 2 ). (2.8) 
Now the result is an immediate consequence of the inverse function theorem. □ 

3 Asymptotic stability and completeness 

This section is devoted to a proof of Theorem ll.il We will first estimate the nonlinearity 
in Subsection 3.1, and then prove the theorem in the subsequent subsection. 



3.1 Estimates on the nonlinearity 

Before starting the proof of our theorems, we establish some nonlinear estimates, first for 
the pointwise nonlinearity, and then for the convolution nonlinearity. 

(I) Pointwise nonlinearity 

Our assumption (|1.25l) implies that for k = 0, 1, 2 

k 

\D k g(z)\ < \g {3) (\z\)z^ k \ < \z\ 7 ^ k + \zf- k . (3.1) 

3=0 

The nonlinear term F%, defined in ()1.15j) . can be expanded by the mean value theorem as 

F 2 (v) = 9(Q + V)~ g(Q) - d° £ g(Q + er,) = f (1 - e)d 2 £ g(Q + e V )de. (3.2) 

Jo 

Then we can estimate it as 

|F 2 |< sup |^(Q + e»7)|H 2 <(l + |Q| + H) 4 H 2 . 

0<e<l 

ll^2|| L i + icc < (l + ||<9||L6 + ||r/|| L 6) 4 N|i 6 . (3.3) 
We also need to estimate g(Q + 77) — g(Q). By using the generalized Holder inequality, 

\\g(Q + v)~ 9iQ)\\ w w S IKAKQ + v)- d 9 (Q))vq\\ l6/5 + \\d 9 {q + r,)v v \\ L e /6 
< (WQWv + IMIl 2 ) 1/3 IMM|vq|| l2 + (HQlUe + ||»7IU«) 3 II»/IIl«I|vq|U8 
+ (\\Q\\% 3 + H\% 3 + HQIIi. + WvWfoWVvWifl 

^ C(\\Q\\hz + \\v\\m)\\v\\w^, 
where C(s) < s 4 / 3 + s 4 . 

(II) Convolution nonlinearity 

The nonlinear term F 2 has the following form: 

F 2 = Q$ * |??| 2 + 77$ * (2Re(Qrj) + |r/| 2 ). 
By the generalized Young inequality in Lorentz spaces, we have under the assumption ()1.26j) . 

ll^|| L i +L oc < (IIQILe + l|r?ILe) IMUm . (3.5) 

This is the only place where we need the Lorentz space 

As for g(Q + 77) — g(Q), its gradient V(g(Q + rj) — g(Q)) is expanded into a sum of 
trilinear forms where one of three functions has the derivative and at least one of them is 
77. By the generalized Young inequality, we have 

||$ * (^1^2)^31^6/5 ^ Plln+LS/a.ooH^i) ||i2 ||V> CT (2) IUviL 6 IIVV(3)IU6, (3.6) 

for any permutation a. So we may put an 77 or Vr) in L 6 , another function without derivative 
in I? n L 6 , and the remaining one in I? . Hence we obtain 

\\g(Q + v)-g(Q) \\ w ^ < \\v\\w^(\\v\\m + \\Q\\m) 2 - (3.7) 
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3.2 Asymptotic stability and completeness 

Now we prove our main result, Theorem ll.il 

Let tp(x,t) solve the nonlinear Schrodinger equation (|1.1|) with initial data 

V>(0, •) = ipo, ||V>o|Ihi( K 3) < l. 

It is easy to prove local well-posedness in H 1 by using the Strichartz estimate (|1.22|) (the 
discrete spectral part does not bother us on finite intervals). The unique solution thereby 
obtained belongs to LfH 1 n L^W 1 ' 6 . 

Our argument below will yield time-global a priori estimates, so that the solution tp 
exists and remains small in H 1 for all time. More precisely, we take 5' > much smaller 
than any 5 in the previous lemmas, and take the initial data ipo such that 

<5'<6. (3.8) 

We will show that 

\\i>\\L?>m[o,T] < $ implies that \\ip\\ L °°m[o,T] < <V 2 > ( 3 - 9 ) 

for any T > 0, provided 5 and 5' were chosen sufficiently small. Then, by continuity in 
time, this bound and the solution together extend globally in time. In the argument below 
we will not explicitly specify the time interval. The assumption HV'Hl^hi < & allows us to 
use all of the previous lemmas. 

By Lemma 12.31 we have the decomposition 

iP = Q[z{t)]+r]{t), rieOQz]. 

The equation for r/ is given in (|1.13j) . We now derive the evolution equation for z(t). 
Differentiating the relation (irj, DjQ[z]} = with respect to t and plugging equation (|1.13JI 
into that, we obtain 

= (Hri + EQ- iDQz + F 2 ,DjQ) + (irj, DjDQz), 

where H, E and Q all depend on z (but this dependence is dropped from the notation). By 
the symmetry of H and (|1.19|) . we have 

(Hrj,DjQ) = (rj, HDjQ) = (r],Dj(EQ)) = (77, EDjQ) = (ir/, EDjDQiz), (3.10) 

where we used {ir], DQ) = 0. By (|1.11|) . we have 

(EQ - iDQz, DjQ) = (DQ(iEz + z),iDjQ). (3.11) 

Thus we obtain 

{(iDjQ,D k Q) + (i V ,D j D k Q))(z + iEz) k = -(F 2 ,D J Q). (3.12) 

fe=l,2 
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The matrix on the left hand side is the Jacobian matrix in (|2.7|) . and so is estimated as 

(iDjQ, D k Q) + (i V , DjD k Q) =j-k + o{8). (3.13) 
Inverting this matrix, we obtain 

\z + iEz\ < \(F 2 ,DQ[z})\ < ||F 2 || L1+LO o (3.14) 
at any t. Applying the estimates ()3.3|) and ()3.5j) . we obtain 

\\Z + IEZ\\ L 2 < ||??|| L oo L 6,2||r/|| L 2 L 6,2 < Wnh^m 11^11^^6,2, (3.15) 

and ||i + iEz\\L°° < IMl!°°.ffi> where we used the Sobolev embedding H 1 C L 6 ' 2 . 
Next, we estimate rj by writing the equation (|1.13|) in the form 

id t r] = (-A + V)rj + F (3.16) 

with 

F := g{Q + 77) - g{Q) - iDQ{z + iEz). (3.17) 

Denote r\ c := P c r\ where P c = 1 — 4>o((j)o, •) is defined in Q1.21JI . The Strichartz estimates 
applied to (|3.16|) and Lemma 12.21 yield 

\\v\\x < \\Vc\\x < ||»7(0)||tfi + \\P C F\\ L 2 W i, 6 /5 < \\t/Jo\\m + \\F\\ L 2 W i, 6 /B, (3.18) 

where 

X := LfH 1 n L 2 VF 1>6 n L 2 L 6 ' 2 . (3.19) 
By Lemma 12. II and the estimates (|3.4|) and (j3.7|) in the previous subsection, we obtain 

\\ F \\L 2 t w^ 6 /s ~ \\z + iEz\\ L 2 + 6\\r}\\ L 2 W x,6. (3.20) 

From (f3~T5|) . (|3~TH|) and (|3~27I1) . we deduce that 

||77||x + p + *£*||j£ 2 + llFllif^i.e/e < HVoki < (3.21) 

if we take 5 sufficiently small. Choosing 5' even smaller, we obtain the desired bootstrapping 
estimate (|3,9|) . and so the solution, as well as all the estimates, extends globally. 
Moreover, we have 

\\d t \z\\\ L i < \\z + iEz\\ L i < \\vf L 2 L e,2 < Ho\\m, (3.22) 

so \z(t)\ and E[z(t)] = E[\z(t)\] converge as t — > 00. 

Finally, we prove that 77 is asymptotically free. We have the integral equation 



Vc (t) = e^ A ^) 



c (0) -i I e- is ^- v ^P c F{s)ds 
Jo 
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(3.23) 



By the Strichartz estimate, for any T > S > we have 

r-T 



I e - is (A-V)p cF ^ di 

Js 



< \\F\ 



H 1 



L 2 W 1 - 6 / 5 [S,T] 



0, 



as T > S — > oo, by the Lebesgue dominated convergence theorem, and the finiteness of 
II-^IIl 2 W 1 ' 8 / 5 (o oo)- Thus the integral in (|3.23j) converges in H 1 as t — ► oo, and we obtain 

;>oo 

lim e- itiA - y) Vc(t) = rfc(O) - i / e- ls{A ~ v) P c F(s)ds =: r/+. 

Jo 

In particular, rj c (t) converges to weakly in H . Then Lemma 12.21 implies that T]d(t) = 
(R[z{t)] — I)rj c {t) converges to strongly in H 1 . Therefore we conclude that 



rj(t) - e^-^TH 



H 1 



0. 



□ 



4 Nonlinear wave operator 

In this section we prove Theorem 11.21 We will construct the desired solution by first 
assigning the asymptotic data at large finite time T and then taking the weak limit as 
T — > oo. Recall that moo, 11??+ II// 1 — For an Y T > 0, we define ifj T to be the solution of 
(|1.1|) with the initial condition 

V{T) = Q[m 00 ] + ^ A ^r ]+ . 

Theorem II .11 implies that tp T is globally defined uniquely and satisfies 



\L 2 t W 1 > fi nL°°H 



\dt\z 



lL!(T,oo) 



<Uz + iEz) T \\ L1{Tj0o) <5 2 



(4.1) 



where the solution can be decomposed according to Lemma 12. HI as 

^ = Q[z T ] + r, T , V T GM C [/]. 
Let £t '■= e~ lT ( A ~ v ^ P C 7] T (T) . Then we have the integral equation for any T > 

P c r, T (t) = e^ A - y kT ~ i f e l ^ A -^P c F T (s)d S , 

where F is as defined in (ETlTJ) . By the same argument as in Subsection 13.21 we deduce 
that for any S <T, 



^WlIW^S.oo] ~ W^V 1 \\ L 2 W 



'• e [S,oo] ~ 



F 



LjW>- 6 / 5 [S,oo] ~ IK Z + lEz Y \\l 2 [S,oo] + ^ H r ? T |ll,fW 1 > 6 [,S,oo] ~ * 11^ IlLfWM^oo] 



lIvv^s,™] 



+ \\F 



| L 2 W 1,6/S [5 j00 ] 



(4.2) 
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Therefore, when 5 is sufficiently small, we have 

MA-V)£ T 



\n T \\ < 



Applying the Strichartz estimate once again, we get 



Per] 1 



MA-v) 



< 6 



(4.3) 



Now we take the limit T — ► oo. First we check the convergence of the data at t = T. 
Denote zt '■= z T (T) and tjt '■= tl T (T)- Since zr is bounded, it converges to some 6 C 
along some subsequence. Then we have 

= {irfr,DQ[z T ]i(moo - z^)) 



= {i(Q[moo] - QlzooDyDQlzool^moo - z^)) = |moo - Zoo| 2 (l - 0(6)), 

where the second equality is by taking limits, and we used Lemma 12. II that rjx = Q[m 
Q[zt] — e lT ^~ v ^r] + and that e lT ^~ v ^r] + converges to weakly in H 1 . Thus 

z T (T) 



(4.4) 



m. 



£t -> V+ in H 1 , 



(4.5) 



as T — > oo (without restriction to a subsequence). 

Now we proceed to convergence for all time. By 1)4.1(1 . z T is equicontinuous on R, and 
so is \z T \ on the extended real line [—00,00]. From the equations, rj T is equicontinuous 
in Ci^w-H 1 ), so is P c r] T in C([-oo, 00], w-H 1 ) by fOJ). Then by Lemma O r] T = 
R[z T ]P c r] T is also equicontinuous in C([— 00, 00], w-H 1 ). Therefore rj T and z T are convergent 
along some subsequence in the following topology: 

- T - ~°° in (C° n L°°)(w-F 1 ) n w-L?^ 1 ' 6 , 



»7 



^0 



z J -^z°°, inC 
This implies the convergence of i(j t itself: 



in L c 



= Q[z T ] + r] T -> Q[z°°] + 77°° =: V c 



(4.6) 



(4.7) 

in C(w-.ff 1 ) fl w-L^ W ' on any finite time interval. Extracting a subsequence if necessary, 
we may assume that the nonlinearity g(ip T ) also converges in w-L^VF 1 ' 6 / 5 on any finite 
time interval. Then the local convergence of tp T in C(L P ) for p < 6 implies that the limit 
of the nonlinearity is the desired g^ 00 ). Hence we deduce that ip°° is a solution to 1(1.1(1 
belonging to C(R; -ff 1 ) n L^ oc (lF 1,6 ). From the uniform convergence of \z T \ to \z°°\, and the 
convergence of z T (T) we have 



lim |z°°(T)| 



lim |z T (T)| = m c 

T — >oo 



From convergence of £r, 1(4.3(1 and the weak convergence uniform in time, we get 
||iV°(t)||#i ->0, ast^oo 



Pert 



< 



H/ 1.6 [5,oo] 



0, as S 



00. 



(4.8) 



Thus ip°° is a solution with the desired asymptotic profile. □ 
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5 Examples of slow decay of dispersion 

In this section we prove Theorem 11.31 

For a fixed ball BcK 3 , choose £o G H 1 satisfying 

UoWm = 1, II£o||l 2 (b) > °> (0o,6) = Oo,£o)l 2 (£0 = 0, (5-1) 

where the first inner product is in L 2 (M 3 ) as before. The constants below depend on B and 
£o- We use a small parameter e > to control the size of solution. 

We define an increasing sequence of times Tj inductively as follows. Let Tj > 1 and for 
j > 1, assume that we have defined 2*. for k < j. By the Strichartz estimate, there exists 
T > maxfc<j Tfc such that 

£ ll^ (A - y)(i - Tfc) eoL^x,a [T)0o] < £ 2^ (5.2) 
fe<i, ± 

and 

sup/(t) < e 2 2^'. (5.3) 
Then we can choose Tj > T such that 

£ ||^ A - v )P>-^&|| W rx..<e2^. (5.4) 
fc<j, ± 

We define the final data by 

V+ :=£ £ 2-V(^)H^ , 
i>o 

and the asymptotic profile of the dispersive part is given by 

m :=J2&~ j j {A ~ V){t ~ Tj ko. (5.5) 

3 

Let ip(t) = Q[z(t)] +i](t) be the solution furnished by Theorem 11.21 corresponding to 77+ 
and moo £ [0, e]. By (|4.3|) and (|5.2|) we have 

\\p c v(Tj) - ru{Tj)\\w- Z % ( 5 - 6 ) 

By (|5.4[) we have 

||r ?l (r j )- £ 2-^ || i6 <e 2 2->\ (5.7) 

Then by Lemma 12.21 we have 

\\PdV(Tj)\\ L e < e\\P c r,(Tj)\\ L e < e 2 2-K (5.8) 
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Thus we obtain 

UTj) - e2-^o\\mB) < WviTj) - e2~^ Q \\ L6 < e 2 2^ . (5.9) 
By Lemma 12. II we have 

IMzft)] - Q[z] - (z(Tj) - z)Mi?(B) < oieMTj) - z\. 
Since (fro > everywhere, we may assume, by choosing e sufficiently small, that 

\\Q[z(T 3 )\ - Q[z] - (ziTj) - z)fo\\v>(B) < eWW) - z)<h\\v>{By (5-10) 
Since an£ i 00 are orthogonal in L 2 (B), we obtain from (|5.9|) and (|5.10|) . 

WvW + QMW-QWWvw >e2-3U \\ LHB) -e 2 2-3 > e2^, 
provided e is sufficiently small. Thus we obtain by (|5.3() 

inf \\ip(Tj) - Q[z]\\ L2{B) > e2'i > VfiTj). 

□ 

6 Appendix: nonlinear bound states 

In this appendix we prove Lemma \2. II 

For the linear potential V, we may weaken Assumption 1 to include only those parts 
which are relevant for existence of nonlinear bound states: 

Assumption 1': We suppose V G L 2 + L°° with || V- Wl^+l^ (W > R) ~> as R ~> °°; an d 
that eo < is a simple eigenvalue of — A + V (we do not need the Strichartz estimates, and 
we do not need eo to be the only eigenvalue, or even the ground state). 

We need the nonlinearity g to be just super-quadratic. Thus we may replace Assumption 
2 by the following: 

Assumption 2': g is as in Assumption 2, but in the pointwise case is only required to 
satisfy (instead of l|1.25j) ) 

g"(s) = o(l) as s£K^0 + . 

Proof of Lemma 12.11 under Assumptions 1' and 2': 

For each z, we look for a solution Q = z(fro + q and E = eo + e' of ()1.6|) with ((fro, q) = 
small and e' € R small. Then 

(-A + V)q + g(z(fr + q) = e'(z(fr + q) + e q. 
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Taking projections on the 4>o and 4>q directions, we get 

e'z = (4> ,g(z(/) +q)) (6.1) 
Tq = -P c g{zfa + q)+ e'q, (6.2) 

where we denote T := —A + V — eo- The right sides are of order o(z). We will use a 
contraction mapping argument to solve for q = o(z 2 ) in H 2 and e' = o(z) uniquely, for 
sufficiently small z. Differentiating by z, we obtain the equations for higher derivatives: 

zDe' + e'J= (<f> ,Dg(Q)), 

zD 2 e' + JDe' = (cb ,D 2 g(Q)), 

TDq = -P c Dg(Q) + qDe' + e'Dq, 

TD 2 q = -P c D 2 g(Q) + qD 2 e' + DqDe + e'D 2 q, 

Dg(Q)=g'(Q)(Jfo + Dq), 

D 2 g{Q) = g"(Q)(Jfo + Dqf + g'(Q)D\ 

where we have omitted subscripts for D and J := Dz = (1, i), and some constant coefficients. 

Assumption 1' implies that —A + V is self-adjoint on 1? with domain H 2 , so 0o £ H 2 . 
Furthermore, the assumption ||T/_ ||(l 2 +l° c )({|x|>r}) ~ * implies that — A + y is a relatively 
compact perturbation of —A + V + . So by the Weyl theorem, the essential spectrum of 
—A + V is contained in [0, oo). In particular, eo is an isolated point of the spectrum. Since 
it is a simple eigenvalue, we have 

T- 1 :L\^H\ bounded, (6.4) 

where H\ and L\ denote the Sobolev spaces restricted to the orthogonal complement of 

00- 

Now we can solve the equations in the closed convex set 

K := {(q, e') E Hi x R | \\q\\ H 2 < \z\ 2 , \e'\ < \z\} (6.5) 

for sufficiently small z £ C (the case z = is trivial). Indeed, we define the map M : 
(9o, eo) i-> (qi,e[) by 

90 ■= g(z(/)o + qo), 

zd x := (4> ,go), (6.6) 
qi :=T-\-P c g + e' o qo). 

Then under Assumption 2' we have the easy estimates 

Wi\ < Hflbllff-a < o(\\z(f> + q Q \\ 2 H 2) < o(z 2 ) 
I|9i||h2 ^ ll5o||i2 + MIMI/p < 0(2; ), 
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which implies that M maps K into K. Let (e'j +2 ,qj+2) '■= M(e?j,qj) and gj := g(z4>o + qj) 
for j = 0, 1. Similarly, we can estimate the difference 

k(4 - 4)1 ~ llfo - 9i\\h-2 < o(\\z4>o + qjWn^ho ~ Qi\\h 2 < o(z)\\q - qi\\ H *, 
\\Q2 ~ 93||h 2 ^ llffo - 9i\\i? + |e' - eiHMlL 2 + killko - QiWm (6-8) 
< IzKleo-eil + Hgo-ftll^) 

which implies that M is a contraction on K, and hence has a unique fixed point in K. 

Suppose now there is a solution (q,e') in the class K' = {(q,e') : \\q\\ H 2 < 7, |e'| < 7}. 
We have 

lklli?2 < ||ff(20O +g)|| L 2 + |e'| \\q\\ L 2 < o(l)\z\ 2 +o(l) ||g|| H2 . 

Thus ||^|[^2 < o(l)|z| 2 . It follows that \ze'\ < ||<7(z0o + <7)||#-2 < o(z 2 ) and hence (g,e') G 
i^. This shows the uniqueness in the class if'. 

Let (q,e') be the unique solution and Q := Z(/>o + q. Since the equation becomes real- 
valued when z£l, the unique solution is also real- valued. 

By the same argument as above, we have 

\z\\De'\<o(z) + \\Dg(Q)\\ H - 2 , 

\\Dq\\m < \\Dg{Q)\\l* + \z\ 2 \De'\ + \z\\\Dq\\ H 2, (6.9) 
\\Dg(Q)\\ L 2 < (||Q|| H2 )|| J0o + Dq\\jp < o(z)(l + \\Dq\\ H 2), 

which imply that De' = o(l) and ||Z)g||^a = o(z). Similarly we have 

\z\\D 2 e'\ <o{l) + \\D 2 g{Q)\\ H -2, 

\\D 2 q\\ H 2 < \\D 2 g{Q)\\ L 2 + o(z 2 )\D 2 e'\ + o{z)\De'\ + o(z)\\D 2 q\\ H 2, (6.10) 
\\D 2 g(Q)\\v < o(l)\\J<fo + Dq\\ 2 H2 + o(z)\\D 2 q\\ H 2, 

which imply that |-D 2 e'| = o(l/z) and jj-D 2 ^)^ = o(l). 

Next we establish the estimate in W 1 ' 1 . Actually, following [Q, we will obtain an expo- 
nentially weighted energy estimate. Let 

£ftM) := (V^,V0) + J V^cpdx 

denote the bilinear form associated to —A + V. It is defined on H 1 x H . Set 

b := lim inf {£(</>, ip) \ 92 € H 1 , \\(p\\ 2 = 1, = for |x| < R}. (6.11) 

Suppose 6 < 0. Then there exists a sequence (fR satisfying ||v_r||2 = 1, ¥?-r( x ) = f° r 
\x\ < R, and £(</?_/?, ^r) < 5 for some fixed 5 < 0. It is easy to check that ipn is bounded in 
JJ . Since it converges weakly to as R — > 00, by the assumption HV-ll/^^oovrui^u — ► 0, 
the negative part J V- \(fji\ 2 dx of the energy converges to 0, a contradiction. Thus b > 0. 
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In other words, there exists 5(R) with 8(R) — > b > as R — > oo, such that for any <p £ if 1 
satisfying <p(x) = for \x\ < R, we have 



|2 
12- 



£(^)> 

Now we apply this inequality with localized exponential weight xr defined by 

' e e(\x\-R) _ Xj ^ < | x | < 2 _ R ) ) 



e e(aR-|as|) _ j ( 2j R < |x| < 3R), 
else 



with some fixed small e > 0. Then we have 

|Vxr| < e( X R + l). 
From (|6.12|) . we have for any tp G H 1 , 

8(R)\\XR<P\\1 < £{XRV,XRV) = ZiXRViV) + J \^XR\ 2 dx 

= {XRV,T<p) + eo\\XR'P\\l + j \^XR\ 2 dx. 
By using Q6.14JI . we obtain for sufficiently large R and small e, 

\\XRV\\l ^ (XR<P,Tip) + £ 2 ||^ll2- 
The Sobolev norm also is estimated by 

llx^v^m < ||V( X ^^)lll + ll^vx^lll 

^ £(XR<P,XR<p) + \\XR<P\\l + e2 |l9 J ll2 
< (x 2 R ip,Tip) + £ 2 yg. 

Thus we obtain the key relation 

WxR^Wm ^ {xr<p,xrT<p) + \\<p\\l, 

for any <p € H 1 , sufficiently large R and small e. Now we substitute each of ip = <f. 
into (|6.18|) and use the equations they satisfy. We find 



(6.12) 



(6.13) 



(6.14) 



(6.15) 



(6.16) 



(6.17) 



(6.18) 
,q,Dq,D 2 q 

(6.19) 



\\XR<t>o\\h < I|0O||| = 1, 

and under Assumption 2' we easily obtain 

WxRqWm ~ WQ^oiQ^Bim^-^iWxRQWmWxRQWm + WQWmWxRQWmWxRMm) + hill 
< (o{z 2 ) + \z\\\xRq\\m)\\xRq\\m + 

(6.20) 



which implies that HX-RSIIh 1 = °( z2 )- Similar estimates hold for XR^q and XR^ 2 q- It 
follows that each of these functions are bounded in W 1 ' . 
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